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Paternal Stress Exposure Alters Sperm MicroRNA Content 
and Reprograms Offspring HPA Stress Axis Regulation 
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Neuropsychiatric disease frequently presents with an underlying hyporeactivity or hyperreactivity of the HPA stress axis, suggesting an excep¬ 
tional vulnerability of this circuitry to external perturbations. Parental lifetime exposures to environmental challenges are associated with 
increased offspring neuropsychiatric disease risk, and likely contribute to stress dysregulation. While maternal influences have been extensively 
examined, much less is known regarding the specific role of paternal factors. To investigate the potential mechanisms by which paternal stress 
may contribute to offspring hypothalamic-pituitary-adrenal (HPA) axis dysregulation, we exposed mice to 6 weeks of chronic stress before 
breeding. As epidemiological studies support variation in paternal germ cell susceptibility to reprogramming across the lifespan, male stress 
exposure occurred either throughout puberty or in adulthood. Remarkably, offspring of sires from both paternal stress groups displayed 
significantly reduced HPA stress axis responsivity. Gene set enrichment analyses in offspring stress regulating brain regions, the paraventricular 
nucleus (PVN) and the bed nucleus of stria terminalis, revealed global pattern changes in transcription suggestive of epigenetic reprogramming 
and consistent with altered offspring stress responsivity, including increased expression of glucocorticoid-responsive genes in the PVN. In 
examining potential epigenetic mechanisms of germ cell transmission, we found robust changes in sperm microRNA (miR) content, where nine 
specific miRs were significantly increased in both paternal stress groups. Overall, these results demonstrate that paternal experience across the 
lifespan can induce germ cell epigenetic reprogramming and impact offspring HPA stress axis regulation, and may therefore offer novel insight 
into factors influencing neuropsychiatric disease risk. 


Introduction 

Stress dysregulation is a common underlying feature in neurode- 
velopmental and affective disorders, with hyperreactivity and hy¬ 
poreactivity of the HPA stress axis reported in schizophrenia, 
autism, and depression (Arborelius et al., 1999; Moghaddam, 
2002; Nestler et ah, 2002; Walker et ah, 2008; Corbett et ah, 2009). 
Parental exposures to stress, infection, malnutrition, and ad¬ 
vanced age have been linked with increased offspring presenta¬ 
tion of neuropsychiatric diseases (Kinney et ah, 2008; Rosenfield 
et ah, 2010; Brown, 2011; Hultman et ah, 2011; Patterson, 2011). 
The HPA stress axis is remarkably vulnerable to environmen¬ 
tal perturbations, and parental experiences likely influence 
disease risk by altering stress pathway regulation (Matthews, 
2002). While maternal insults experienced during pregnancy 
can directly impact fetal development, the mechanisms by 
which lifelong maternal or paternal exposures can affect off¬ 
spring neurodevelopment are less well understood (Bale et ah, 
2010 ). 
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Epigenetic marks in germ cells are a proposed mechanism by 
which parental environmental exposures shape offspring neuro¬ 
development (Jirtle and Skinner, 2007; Guerrero-Bosagna and 
Skinner, 2012). When examining mechanisms of transgenera- 
tional transmission, studies of maternal experience must con¬ 
sider the complex maternal-fetal/neonatal interaction, including 
changes in the intrauterine environment, lactation, and maternal 
behaviors (Champagne, 2008). However, paternal studies in ro¬ 
dent models, where the sire does not participate in rearing of the 
offspring, allow for relative exclusion of factors outside of the 
germ cell contribution, providing a homogenous and easily ac¬ 
cessible tissue for examining epigenetic marks (Rando, 2012). 

Research from rodent models of male exposure to dietary 
challenge, drugs of abuse, or social defeat support the transmis¬ 
sion of paternal experiences to offspring through epigenetic 
marks in sperm (Carone et al., 2010; Guerrero-Bosagna et al., 
2010; Ng et al., 2010; Dietz et al., 2011; Vassoler et al., 2013). 
Although mature sperm lack the machinery necessary for chro¬ 
matin remodeling, spermatozoa RNA populations, histone mod¬ 
ifications, and methylation patterns reflect heritable responses to 
environmental challenges that may occur during spermatogene¬ 
sis (Godmann et al., 2009; Skinner and Guerrero-Bosagna, 2009). 
Of specific interest is the sperm microRNA (miR) content, as 
miRs are a mode of dynamic epigenetic regulation that rapidly 
respond to environmental perturbations and likely regulate gene 
expression postfertilization (Giraldez et al., 2006; Rassoulzade- 
gan et al., 2006; Grandjean et al, 2009; Leung and Sharp, 2010). 
Epidemiological evidence from Swedish famine years suggests 
that germ cell susceptibility to reprogramming, and its subse- 
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quent effects on future generations’ disease risk, is greatest during 
the slow growth period of preadolescence (Kaati et al., 2002, 
2007; Pembrey et al., 2006). 

Therefore, to examine the effect of paternal stress exposure 
across the lifespan on offspring HPA axis regulation, male mice 
were exposed to 6 weeks of chronic stress, before breeding, either 
throughout puberty or only in adulthood. Offspring stress reac¬ 
tivity and gene array analysis of offspring stress regulatory brain 
regions, the the paraventricular nucleus (PVN) and bed nucleus 
of stria terminalis (BNST), were subsequently examined. Pater¬ 
nal sperm miR content was assessed as an indication of potential 
epigenetic marks associated with the reprogramming of offspring 
physiology or behavior. 

Materials and Methods 

Animals 

All mice bred for paternal stress studies were virgin, in-house C57BL/6: 
129 FI hybrid mice. The hybrid vigor of this background strain provides 
a reproducible balance of stress responsivity, behavioral performance, 
and maternal care (Mueller and Bale, 2006,2007,2008; Morgan and Bale, 
2011). All mice were housed in a 12 h light/dark cycle with ambient 
temperature 22°C and relative humidity of 42%. Food (Purina Rodent 
Chow; 28.1% protein, 59.8% carbohydrate, 12.1% fat) and water were 
provided ad libitum. All studies were performed according to experimen¬ 
tal protocols approved by the University of Pennsylvania Institutional 
Animal Care and Use Committee, and all procedures were conducted 
in accordance with the NIH Guide for the Care and Use of Laboratory 
Animals. 

Chronic variable stress 

Administration of chronic variable stress was performed as described 
previously (Mueller and Bale, 2006). Briefly, males were randomly as¬ 
signed to a control nonstressed group (n = 13) or to one of two experi¬ 
mental groups, receiving 42 d of chronic variable stress throughout 
puberty (PN28-70; n = 13) or in adulthood (PN56-98; n = 13). Seven 
different stressors were randomized and administered one per day. The 
order of seven stressors varied across weeks. The duration of stress (42 d) 
was selected to encompass a complete round of spermatogenesis, expos¬ 
ing sperm to stress at all stages of maturation (Oakberg, 1956). Stressors, 
selected because they are nonhabituating, do not induce pain, and do not 
affect food or water intake, included the following: 36 h constant light, 15 
min exposure to fox odor (1:5000, 2,4,5-trimethylthiazole; Acros Organ¬ 
ics), novel object (marbles) overnight, 15 min restraint in a 50 ml conical 
tube, multiple cage changes, novel 100 dB white noise (Sleep Machine; 
Brookstone) overnight, and saturated bedding overnight. 

Breeding scheme 

Following completion of chronic variable stress, all males were housed 
with a single female for 14 d to remove the acute effects of stress on male 
behavior and mature sperm stored in the epididymis. Males were then 
bred with naive females for 1-3 nights. To minimize male-female inter¬ 
actions, which may influence maternal investment or care and subse¬ 
quently affect offspring development (Curley et al., 2011), pregnancy was 
established by confirmation of a copulation plug (identified each morn¬ 
ing within 1 h after lights on) and signaled immediate removal of the 
female to her own cage containing a nestlet. Breeding resulted in litters 
sired from control (n = 8), pubertally stressed (n = 10), and adult 
stressed (n = 9) males. 

Characterization of offspring physiological and 
behavioral phenotypes 

Same-sex littermates were split into three groups for physiological and 
behavioral characterization (maximum of one male and female per litter 
per group). Groups were examined for either FfPA stress axis responsiv¬ 
ity and prepulse inhibition (PPI), performance on the tail suspension test 
(TST) and Barnes maze, or performance in the light-dark (LD) box and 
postmortem molecular measures. A minimum 10 d recovery period sep¬ 
arated tests. 


HPA axis responsivity. Plasma corticosterone was measured following 
an acute 15 min restraint in a 50 ml conical tube. Testing occurred 2-5 h 
after lights on. Tail blood from adult mice (n = 6-9 per group) was 
collected at onset and completion of restraint (0 and 15 min, respec¬ 
tively) and 15 and 75 min after the end of restraint (30 and 90 min, 
respectively). Tail blood collection requires <30 s to complete. Samples 
were immediately mixed with 50 mM EDTA and centrifuged 10 min at 
5000 rpm. Plasma was collected and stored at — 80°C until analysis. Cor¬ 
ticosterone levels were determined by 125 I-corticosterone radioimmu¬ 
noassay (MP Biomedicals) following kit instructions. The minimum 
detection limit of the assay was 7.7 ng/ml, and the intraassay coefficient of 
variation was 7.1%. 

Preptdse inhibition. PPI of acoustic startle was recorded in SR-LAB 
startle chambers (San Diego Instruments) as described previously (Geyer 
and Dulawa, 2003). All test sessions were conducted 2- 6 h after lights off. 
A PPI test session began with 5 min acclimation to background noise (65 
dB) followed by five consecutive acute habituation tones (40 ms dura¬ 
tion; 120 dB). Ten repetitions of each of the following trial types were 
then presented in pseudorandomized blocks: startle pulse only (40 ms; 
120 dB), no stimulus (65 dB; background), and prepulses +4, +8, or 
+16 dB above background (20 ms; 69, 73, or 81 dB) preceding the startle 
pulse (40 ms; 120 dB) by 100 ms. Intertrial intervals averaged 15 s. Acous¬ 
tic startle response was defined as the peak startle magnitude recorded 
during 65 consecutive 1 ms readings following the startle pulse onset. The 
percentage PPI for each prepulse intensity was calculated for each animal 
(n = 5-6 per group) as [1 — (average response to prepulse + startle)/ 
(average response to startle only)] X 100 before group averages were 
compared. Animals with average response in startle only trials of <50 
mV were excluded from analysis. 

Tail suspension test. The TST was performed as described previously 
(Steruet al., 1985). Briefly, adult mice (n = 7-10 per group) were secured 
to a rod by adhesive tape placed 1 cm from the tip of the tail and sus¬ 
pended 50 cm above the bench top. Testing occurred 2-5 h after lights 
on. Total immobility time, defined as the absence of movement except 
whisker twitches and respiration, was recorded over the 6 min test ses¬ 
sion. Testing was recorded digitally and analyzed manually using ANY- 
maze v4.75 software (Stoelting). 

Barnes maze. A modified Barnes maze was performed as described 
previously (Mueller and Bale, 2007). Briefly, adult mice (n = 7-10 per 
group) were trained on two trials per day for 3 d (4 h separation between 
trials) to find a target escape box. The location of the target escape box 
under 1 of 24 holes evenly spaced around the perimeter of a 90 cm 
circular disk remained constant throughout training. Latency to identify 
the target escape box by nose poke or head deflection was determined for 
each trial. Each trial ended when the mouse identified and entered the 
target box or after 4 min of elapsed time. If the target box was not 
successfully entered within 4 min, the investigator guided the mouse into 
the target box, and a latency of 240 s was assigned. Testing was digitally 
recorded and analyzed manually using ANY-maze v4.75 software 
(Stoelting). 

Light-dark box. The light-dark box test was performed as described 
previously (Bale et al., 2000; McEuen et al., 2009). Adult mice (n = 6-8 
per group) were placed in the light compartment at the beginning of the 
10 min test session. Light intensities were set at 5 lux in the dark com¬ 
partment and 300 lux in the light compartment. All testing occurred 2-5 
h after lights off. Total time spent in the light compartment and the 
number of light-to-dark transitions were analyzed with ANY-maze v4.75 
software (Stoelting). 

Mechanistic evaluation of blunted stress responsivity 
Acute selective serotonin reuptake inhibitor. Activation of the HPA axis by 
citalopram, in the absence of any additional stressor was examined as 
described previously (Goel et al., 2011). Plasma corticosterone was mea¬ 
sured in offspring (n = 5-8 per group) as described above, except that 
animals were injected with 15 mg/kg citalopram (Sigma-Aldrich) or ve¬ 
hicle (0.9% saline) at time 0 min and were not restrained. 

Adult tissue collection. Adult male and female offspring of paternal 
stress or control sires (n = 6-8 per group) were rapidly decapitated 
under isoflurane anesthesia at 22-24 weeks of age. Whole brains were 


Rodgers et al. • Transgenerational Impact of Paternal Stress 


J. Neurosci., May 22,2013 • 33(21):9003-9012 • 9005 


removed, immediately frozen on dry ice, and stored at — 80°C. The pitu¬ 
itary gland and left adrenal gland were removed, flash frozen in liquid 
nitrogen, and stored at — 80°C. Approximately 5 X 10 6 sperm were 
isolated from the caudal epididymis of adult males into PBS with 1% BSA 
at 37°C through a swim-up assay. After settling for 30 min, sperm- 
containing supernatant was centrifuged for 5 min at 4000 rpm. Sperm 
pellets were stored — 80°C. 

Peripheral gene expression. Whole pituitary and left adrenal glands 
(n = 6-8 per group) were bath sonicated for 7 min (30 s on/off cycles) 
in TRIzol reagent (Invitrogen), and total RNA was isolated according 
to manufacturer’s protocol. cDNA was transcribed using the High- 
Capacity cDNA reverse transcriptase kit (Applied Biosystems). Changes in 
pituitary corticotrophin-releasing factor receptor 1 (CRFrl; NM_007762.4, 
Mm00432670_ml), pituitary proopiomelanocortin (POMC; NM_ 
008895.3, Mm00435874_ml), adrenal melanocortin 2 receptor (Mc2r; 
NM_008560.2, Mm00434865_sl), and adrenal 1 l-j8 hydroxysteriod dehydro¬ 
genase 1 (11J3HSD-1; NM_001044751.1, Mm00476182_ml) were measured 
by quantitative real-time PCR using TaqMan gene expression assays (Ap¬ 
plied Biosystems). Samples were run in triplicate for the target genes and 
endogenous controls (pituitary, Gapdh, NM_0080842.2, Mm99999915_gl; 
adrenal, Actb, NM_007393.3, Mm00607939_sl) on the same 96-well plate. 
Analysis was performed by the comparative Ct method, and expression 
levels were normalized to control male offspring (Schmittgen and Livak, 
2008). 

Brain gene expression. Whole brains were cryosectioned at — 20°C. 
Using a hollow needle (Ted Pella), brain regions were micropunched 
according to the atlas of Paxinos and Franklin (2001) as follows: BNST, 
0.75 mm bilateral punches from one 300 pm slice +0.25 to —0.05 rela¬ 
tive to bregma (atlas Figures 25-32); PVN, 1.00 mm punch along the 
midline from two successive 300 pm slices —0.50 to —0.80 and —0.80 to 
—1.10 relative to bregma (atlas Figures 36-40). Punches were bath son¬ 
icated for 2.5 min (30 s on/off cycles) in TRIzol reagent (Invitrogen). 
RNA was isolated by RNeasy kit (Qiagen) and suspended in RNase- 
free water. Total RNA (n = 5-6 per group) was sent to the University 
of Pennsylvania Path BioResource Molecular Profiling Core for Af- 
fymetrix GeneChip Mouse Gene 1.0 ST analysis. Gene set enrichment 
analysis (GSEAv2.0.7; Broad Institute, Cambridge, MA) was used to 
interpret patterns of gene expression (Mootha et al., 2003; Subrama- 
nian et al., 2005). Collections of c2 (curated: Kegg, Biocarta, and 
Reactome), c3 (motif), and c5 [gene ontology (GO)] annotated gene 
sets were obtained from the Molecular Signature Database 
(MSigDBv3.0.1; Broad Institute) available for use with GSEA 
software. 

Analysis of sperm miR environment 

Total RNA (n = 4 per group) was extracted from sperm pellets using 
TRIzol reagent (Invitrogen). Total RNA (100 ng) was reverse tran¬ 
scribed to cDNA using Megaplex RT pool A and B primers and Mul¬ 
tiscribe reverse transcriptase, and then quantified using the TaqMan 
Array MicroRNA card A and B Set v3.0 (Applied Biosystems). Anal¬ 
ysis was performed using the comparative Ct method (Schmittgen 
and Livak, 2008). For each sample, the average of the Ct values of 
mammu6 and sno202 was used as an endogenous loading control. 
Expression levels of each sample were normalized to the average ex¬ 
pression level of male controls. miRs with expression in <50% of 
samples or Ct values >35 were excluded, such that analysis of signif¬ 
icant differences in expression of specific miRs was conducted with 
n = 2-4 per group. The miRWalk database, which identifies putative 
miR target sequences in mRNA 3' -UTR commonly predicted by miR¬ 
Walk, miRDB.org, miRanda, and TargetScan algorithms, was used to 
determine predicted targets of significantly changed miRs (Dweep et 
al„ 2011). 

Statistics 

An investigator blind to animal treatment group conducted all studies 
and analyses. Total corticosterone was analyzed by multivariate ANOVA 
(sex by paternal stress for restraint experiment; sex by paternal stress by 
drug for citalopram experiment) with time as a repeated measure. PPI 



Time (min) 



Time (min) 

Figure 1. Paternal stress experienced throughout puberty or in adulthood elicited stress 
dysregulation in offspring. A, 8, Both male (/I) and female (8) offspring of sires that had been 
exposed to chronic stress throughout puberty or only in adulthood produced less corticosterone 
relative to offspring of control sires following a 15 min restraint (shaded column). Total AUC of 
corticosterone production (shown in the inset) showed a significant decrease in both stress 
groups relative to controls. An expected sex difference in corticosterone production was ob¬ 
served, with females showing a greater response than males. Data are presented as mean ± 
SEM.n = 6 -9 litters per group. *p < 0.05. 


was analyzed by two-way ANOVA (sex by paternal stress) with prepulse 
intensity as a repeated measure. Baseline acoustic startle, TST immobility 
time, and LD box measures were analyzed by two-way ANOVA (sex by 
paternal stress). Latency to find target on the Barnes maze was analyzed 
by two-way ANOVA (sex by paternal stress) with trial as a repeated 
measure. Adrenal and pituitary gene expression analysis was analyzed by 
two-way ANOVA (sex by paternal stress). Data greater than two SDs 
away from each group mean were considered outliers and excluded from 
analysis. Main effects and interactions were further analyzed with 
Tukey’s HSD post hoc test. Significant differences were identified at p < 
0.05. All statistics above were performed using JMP8 (SAS) software, and 
all data are reported as mean ± SEM. Microarray data were modeled as 
“control” versus “paternal stress” in accordance with behavioral and 
physiological data. Thresholds for multiple comparisons were set at 
false discovery rate (FDR) £ 0.05. Phenotype permutations (1000) 
were computed in GSEA to determine the FDR, nominal p value, and 
normalized enrichment score (NES) of each gene set. Gene sets with 
FDRs of <0.25, p £ 0.01, and NESs of >1.6 met significance thresh¬ 
olds. For miR analysis, heat maps and hierarchical clustering were 
performed using MuliExperiment Viewer (TM4.org) with Pearson’s 
correlation, and miRs significantly altered in both pubertal and adult 
stress groups were identified by one-way ANOVA with post hoc Stu¬ 
dent’s t test (p < 0.05). 
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Figure 2. Paternal stress does not alter offspring performance in behavioral assessments or HPA response to an SSRI. A, B, 
Neither male (ll) nor female (8) offspring of stressed sires showed altered PPI of the acoustic startle response compared to control 
offspring. C, PPI max startle was also unchanged by paternal stress in both sexes. D, No significant differences in time spent 
immobile on the tail suspension test were observed between control and paternally stressed male or female offspring, f, F, In the 


Results 

Analysis of offspring physiology 
and behavior 

Offspring HPA axis 

To determine whether paternal exposure 
to chronic stress during pubertal or adult 
windows altered offspring stress respon- 
sivity, we examined plasma corticosterone 
levels following a 15 min restraint stress. 
There were significant effects of paternal 
stress (F( 2 j35 ) = 9.3, p = 0.0006) and sex 
(F ( i, 35 ) = 22.68, p < 0.0001) on cortico¬ 
sterone levels that did not interact 
(F( 2 , 35 ) = 0.58, p = 0.56), and a significant 
effect of time (F ( 3j33 ) = 107.4, p < 0.0001) 
(Fig. 1 A,B). Similarly, there were signifi¬ 
cant effects of paternal stress (F (2>3 5 ) = 5.9, 
p = 0.006) and sex (F (I35) = 18.7, p = 
0.0001) on total corticosterone that did 
not interact (F (2j35) = 0.2, p = 0.86). Total 
corticosterone was lower in paternal stress 
conditions compared to controls and in 
males compared to females {p < 0.05) 
(Fig. 1 A,B, inset). The following statisti¬ 
cal outliers in corticosterone levels at each 
of the four time points were removed: one 
pubertal male and one pubertal female at 
time 0, one pubertal female at time 15, 
none at time 30, and one adult male and 
one adult female at time 120. 

Offspring behavior 

To examine the degree to which blunted 
stress axis responsivity was representative 
of global stress sensitivity in offspring of 
stressed sires, we assessed offspring PPI, 
stress-coping behavior in the TST, spatial 
learning in an adverse environment on the 
Barnes maze, and anxiety-like behavior 
in the LD box. PPI did not differ between 
males and females (F (123) = 0.54, 
p = 0.47) or between paternal stress and 
control groups (F (223) = 0.28, p = 0.75), 
though an expected within-subject effect 
of prepulse intensity was observed, with 
higher prepulses eliciting greater prepulse 
inhibition (F (2>22) = 64.1, p < 0.0001) 
(Fig. 2 A,B). Max startle to 120 dB tone 
was also unaffected by sex (F (1>23) = 0.55, 


Barnes maze spatial learning and memory task, no differences 
were detected in the latency to learn the task. 6, H, In the 
light-dark box, no treatment group differences were de¬ 
tected for time spent in the light or transitions between the 
compartments. There was an overall effect of females spend¬ 
ing more time in the light compartment than males. /, J, Fol¬ 
lowing administration of an SSRI, citalopram, at time 0 (black 
arrow), male and female offspring of stressed or control sires 
exhibited similar stress reactivity. Total AUC of corticosterone 
production (shown in the inset) did not vary with sex or pater¬ 
nal stress group. Data are presented as mean ± SEM. n = 5-9 
litters per group. *p < 0.05. 
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Figure 3. Stress axis-related gene expression in offspring pituitary and adrenal glands was not significantly different between 
control (C), pubertal stress (P), or adult stress (A) groups. A-D, Real-time PCR analysis revealed no effect of paternal stress on 
pituitary CRFrl (71), POMC (B), adrenal Mc2r (f), or 1 1 /3HSD-1 (D). Sex differences were observed, with males showing higher 
expression of CRFrl and 11 /3HSD-1 and lower expression of Mc2r and POMC relative to females. Data are presented as mean ± 
SEM./r = 6-8 litters per group. **p < 0.001 ;***p < 0.0001. 


Table 1. Summary of gene sets meeting GSEA enrichment criteria of FDR < 0.25, 
/>< 0.01,and NES > 1.6 



Total 

PVN (enriched) 

BNST (enriched) 

c2 (curated) 

305 

0 (0.0%) 

0 (0.0%) 

c3 (motif) 

776 

16 (2.1%) 

114(14.7%) 

c5 (gene ontology) 

969 

0 (0.0%) 

6 (<1%) 


Table 2. Gene sets enriched in offspring PVN following paternal stress 

Gene set 

Type 

Size 

NES 

P 

FDR 

NR3CI (GR) 

c3 

107 

-1.683 

0.008 

0.235 

CEBPdelta 

c3 

212 

-1.682 

0.001 

0.220 

F0XD3 

c3 

178 

-1.682 

0.006 

0.207 

miR-522 

c3 

144 

-1.678 

0.008 

0.203 

NKX2-5 

c3 

109 

-1.665 

0.008 

0.196 

PAX8 

c3 

88 

-1.661 

0.004 

0.194 

TITF1 

c3 

212 

-1.657 

0.002 

0.194 

PIT1 

c3 

195 

-1.649 

0.010 

0.201 

DDIT3 

c3 

209 

-1.641 

0.004 

0.200 

NFI 

c3 

226 

-1.633 

0.001 

0.193 

F0XI3 

c3 

165 

-1.625 

0.010 

0.193 

IPFI 

c3 

220 

-1.617 

0.001 

0.190 

AR 

c3 

227 

-1.606 

0.002 

0.203 


Members of c3 gene sets share a binding motif for the designated transcription factor or miR. Gene sets defined by 
conserved sequence only and those not meeting the criteria NES > 1.6, p < 0.01, and FDR < 0.25 were omitted. 


p = 0.46) or paternal stress (-F (2i23) = 0.07, p = 0.93) (Fig. 2C). In 
the TST, there was no effect of sex (F (140) = 2.2, p = 0.12) or 
paternal stress (F (240) = 0.01, p = 0.89) on total time spent im¬ 
mobile (Fig. 2D). Analysis of latency to find the target on the 
Barnes maze similarly revealed no effect of sex (F ( 1>29 ) = 0.66, p = 
0.78) or paternal stress (F (229) = 0.41, p = 0.66), and latency to 
find target decreased across trials for all groups (F (5j25 ) = 44.83, 
p < 0.0001) (Fig. 2E,F). In the LD box, females spent signifi¬ 


cantly more time in the light than males 
(-F(i,32) = 5-9, p = 0.02), but there was no 
effect of paternal stress (F (232) = 1.5, p = 
0.21) (Fig. 2G). To determine whether 
Barnes maze and LD box results were af¬ 
fected by locomotor ability, light-to-dark 
transitions were measured in the LD box. 
There was no effect of sex (F (132 ) = 3.44, 
p = 0.07) or paternal stress on the number 
of transitions (F (2 32) = 0.20, p = 0.81) 
(Fig. 2 H). Statistical outliers were removed 
only from analysis of maximum startle and 
prepulse inhibition (one adult female) (Fig. 
2 B, C) and from analysis of total time in light 
(one pubertal male and one control female) 
(Fig. 2G). 

Mechanistic evaluation of blunted 
stress reactivity 

Acute selective serotonin reuptake 
inhibitor treatment 

To investigate the potential role of sero¬ 
tonin in mediating changes in offspring 
stress axis responsivity, we evaluated off¬ 
spring corticosterone levels in response to 
acute selective serotonin reuptake inhibi¬ 
tor (SSRI) administration. Citalopram 
administration elicited similar HPA acti¬ 
vation in all groups, as there was a main 
effect of drug (F ( lj52) = 14.6, p = 0.0004) 
on corticosterone level that did not interact with sex (F (li52) = 
0.003, p = 0.96) or paternal stress (F (252) = 0.88, p = 0.42) (Fig. 
27,/). The following statistical outliers in corticosterone levels at 
each of the four time points were removed: at time 0, one pubertal 
vehicle male, one citalopram control female, one pubertal vehicle 
female, and one adult vehicle female; none at time 15 or time 30; 
and at time 120, 1 pubertal vehicle female and one adult citalo¬ 
pram female. 

Peripheral gene expression analysis 

To determine whether changes in peripheral control of the HPA axis 
explain offspring blunted stress responsivity following paternal 
stress exposure, we evaluated gene expression in offspring pituitary 
and adrenal glands. Neither pubertal nor adult paternal stress signif¬ 
icantly altered the expression of genes involved in the offspring HPA 
axis response, including pituitary CRFrl (F (227) = 1.30, p = 0.29), 
pituitary POMC (F (2>27) = 0.13, p = 0.90), adrenal Mc2r (F (226) = 
2.32, p = 0.12), or adrenal 11/3HSD-1 (F (2>26) = 0.55, p = 0.58). 
Consistent with sex differences in the expression of these genes re¬ 
ported previously (Goel and Bale, 2010), significant effects of sex 
were observed: CRFrl was higher in males (F (I27:) = 21.42, p < 
0.0001), POMC was higher in females (F (1>27) = 14.39, p = 0.001), 
Mc2r was higher in females (F (126) = 24.28, p < 0.0001), and 
11/3HSD-1 was higher in males (F (li26) = 16.92, p = 0.003). There 
were no interactions of paternal stress and sex for CRFrl (F (227) = 
0.24, p = 0.79), POMC (F (2>27) = 0.01, p = 0.98), MC2R (F (2 ’ 26) = 
2.0, p = 0.16), or llj3HSD-i (F (2>26) = 0.05, p = 0.95) (Fig. 3 A-D). 

Brain gene expression analysis 

To examine changes in central regulators of the HPA axis following 
paternal stress, we evaluated patterns of gene expression in the PVN 
and BNST with a microarray. By multiple comparison analysis, no 
single gene met both threshold and significance criteria. By GSEA 
analyses, c3 gene sets in the PVN and c3 and c5 gene sets in the BNST 
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Table 3. Gene sets enriched in offspring BNST following paternal stress 


Gene set 

T yP e 

Size 

NES 

P 

FDR 

Gene set 

Type 

Size 

NES 

P 

FDR 

Chromatin modifications (HAT) 




miR-154 family 






CREBI 

c3 

231 

-1.723 

0.002 

0.027 

miR-381 

c3 

182 

-1.902 

0.002 

0.025 


c3 

188 

-1.721 

0.004 

0.027 

miR-323 

c3 

143 

-1.899 

0.002 

0.021 


c3 

233 

-1.704 

0.006 

0.028 

miR-410 

c3 

79 

-1.808 

0.004 

0.023 

ATF2 

c3 

219 

-1.688 

0.004 

0.031 

miR-494 

c3 

143 

-1.746 

0.006 

0.027 

CEBPalpha 

c3 

215 

-1.664 

0.001 

0.033 

miR-409-3p 

c3 

122 

-1.722 

0.001 

0.028 

AR 

c3 

50 

-1.655 

0.006 

0.034 

miR-382 

c3 

67 

-1.675 

0.010 

0.031 


c3 

227 

-1.621 

0.006 

0.038 

m/R-8 family 






EVII 

c3 

206 

-1.602 

0.002 

0.040 

miR-200A 

c3 

46 

-1.907 

0.004 

0.026 

PCAF 

c3 

213 

-1.613 

0.002 

0.038 

miR-200B,C, miR-429 

c3 

414 

-1.854 

0.002 

0.020 

SREBFI 

c3 

393 

-1.724 

0.001 

0.028 

miR-l4l,miR-200A 

c3 

281 

-1.778 

0.004 

0.025 

STATI 

c3 

213 

-1.724 

0.001 

0.028 

other miR families 







c3 

210 

-1.655 

0.007 

0.034 

miR-l, miR-206 

c3 

260 

-1.827 

0.004 

0.021 

EGRI.EGR2, EGR3 

c3 

240 

-1.764 

0.004 

0.026 

miR-96 

c3 

286 

-1.785 

0.006 

0.026 

ETSI 

c3 

205 

-1.618 

0.010 

0.038 

miR-105 

c3 

154 

-1.783 

0.004 

0.025 

GATAI 

c3 

209 

-1.663 

0.008 

0.033 

miR-IOA,B 

c3 

108 

-1.819 

0.002 

0.021 

MYCN 

c3 

228 

-1.735 

0.002 

0.028 

miR-l28A,B 

c3 

286 

-1.856 

0.002 

0.020 

MYC, MAX 

c3 

215 

-1.728 

0.004 

0.028 

miR-130A,B, miR-301 

c3 

341 

-1.830 

0.004 

0.022 


c3 

212 

-1.690 

0.002 

0.031 

miR-138 

c3 

186 

-1.740 

0.010 

0.027 

MYC 

c3 

225 

-1.738 

0.006 

0.027 

miR-139 

c3 

114 

-1.879 

0.002 

0.020 


c3 

160 

-1.617 

0.004 

0.038 

miR-144 

c3 

175 

-1.746 

0.008 

0.027 

MAX 

c3 

220 

-1.805 

0.001 

0.023 

miR-150 

c3 

79 

-1.713 

0.002 

0.027 

Neuronal survival 






miR-190 

c3 

56 

-1.730 

0.002 

0.028 

E2FI 

c3 

205 

-1.714 

0.004 

0.028 

miR-l99A 

c3 

160 

-1.751 

0.006 

0.026 


c3 

201 

-1.678 

0.006 

0.031 

miR-218 

c3 

343 

-1.788 

0.008 

0.025 

TFDPI 

c3 

205 

-1.686 

0.002 

0.031 

miR-22l,miR-222 

c3 

113 

-1.884 

0.001 

0.022 


c3 

200 

-1.674 

0.004 

0.031 

miR-223 

c3 

85 

-1.779 

0.004 

0.025 


c3 

207 

-1.652 

0.008 

0.034 

miR-23A, miR-23B 

c3 

372 

-1.706 

0.008 

0.028 


c3 

206 

-1.622 

0.008 

0.038 

miR-27A, miR-27B 

c3 

418 

-1.758 

0.001 

0.027 


c3 

186 

-1.658 

0.006 

0.034 

miR-302C 

c3 

219 

-1.926 

0.002 

0.021 

ELKI 

c3 

210 

-1.839 

0.006 

0.021 

miR-330 

c3 

271 

-1.825 

0.004 

0.021 

SRF 

c3 

199 

-1.658 

0.010 

0.034 

miR-335 

c3 

78 

-1.881 

0.001 

0.021 


c3 

190 

-1.654 

0.004 

0.034 

miR-345 

c3 

48 

-1.864 

0.006 

0.020 

Circadian rhythmicity 





miR-361 

c3 

78 

-1.765 

0.010 

0.027 

AHR 

c3 

182 

-1.637 

0.010 

0.035 

miR-362 

c3 

63 

-1.802 

0.002 

0.023 


c3 

126 

-1.617 

0.004 

0.038 

miR-373 

c3 

201 

-1.834 

0.002 

0.022 

ARNT 

c3 

218 

-1.877 

0.001 

0.019 

miR-380-3p 

c3 

88 

-1.901 

0.001 

0.023 


c3 

208 

-1.799 

0.001 

0.023 

miR-412 

c3 

53 

-1.788 

0.006 

0.026 

Other celluar functions 





miR-495 

c3 

213 

-1.810 

0.002 

0.023 

E4FI 

c3 

246 

-1.760 

0.004 

0.027 

miR-498 

c3 

97 

-1.758 

0.010 

0.026 

HNF4A 

c3 

239 

-1.627 

0.001 

0.037 

miR-499 

c3 

65 

-1.784 

0.002 

0.025 

NFE2LI, MAFG 

c3 

255 

-1.621 

0.006 

0.038 

miR-505 

c3 

89 

-1.705 

0.008 

0.028 

PAX3 

c3 

124 

-1.690 

0.008 

0.030 

miR-517 

c3 

43 

-1.893 

0.004 

0.021 

PCBPI 

c3 

221 

-1.640 

0.006 

0.034 

miR-5l9E 

c3 

109 

-1.647 

0.010 

0.034 

POU2FI 

c3 

195 

-1.617 

0.002 

0.038 

miR-520F 

c3 

205 

-1.726 

0.004 

0.028 

SOX9 

c3 

202 

-1.706 

0.001 

0.028 

miR-520G,miR520H 

c3 

195 

-1.720 

0.006 

0.027 

YYI 

c3 

210 

-1.748 

0.001 

0.027 

miR-520A-E, miR-93 

c3 

284 

-1.719 

0.006 

0.027 

ZFPI6I 

c3 

207 

-1.947 

0.002 

0.019 

miR-522 

c3 

144 

-1.824 

0.008 

0.020 







miR-524 

c3 

394 

-1.873 

0.001 

0.019 

Gene Ontology: Transcription 










Helicase Activity 

c5 

43 

-1.762 

0.008 

0.235 

RNA Polymerase II Activity 

c5 

160 

-1.681 

0.001 

0.156 

Nucleoplasm 

c5 

246 

-1.761 

0.004 

0.218 

Isomerase Activity 

c5 

30 

-1.677 

0.010 

0.152 

RNA Helicase Activity 

c5 

22 

-1.748 

0.002 

0.220 

Transcription Factor Binding 

c5 

270 

-1.647 

0.010 

0.157 


Members of c3 gene sets sharea binding motiffor designated transcription factor or miRNA. c3 gene sets sharing the same name represent distinct binding sites for the transcription factor or miR given. Members of c5 gene sets share common 
ontolgy and do not necessarily represent coregulated genes. Gene sets defined by conserved sequence oniy and those not meeting the criteria NES ^ 1.6, p ^ 0.01, and FDR ^ 0.25 were omitted. HAT, Histone acetyltransferase. 
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Figure 4. Paternal stress experienced throughout puberty or in adulthood produced robust changes in sperm miR content. A, Hierarchical clustering of sperm miR expression isolated a single 
dade containing all pubertal and adult stress males, and no controls (bolded red line), n = 4 per group. B, Analyses of significant differences in expression of specific miRs revealed nine that were 
significantly increased in both paternal stress groups. Data are presented as mean ± SEM. *p < 0.05, different from controls; *p < 0.05, different from pubertal stress. 


were enriched following stress (Table 1). Gene sets enriched in off¬ 
spring PVN following paternal stress were sorted according to the 
degree of enrichment (Table 2). Due to the large number of signifi¬ 
cant gene sets in the BNST, enriched gene sets were organized ac¬ 
cording to predicted function or miR family (Table 3). Members of 
c3 gene sets share known or likely regulatory elements in their pro¬ 
moters or 3'-UTRs, and c3 gene sets are named according to the 
transcription factor or miR that binds that conserved sequence. c3 
gene sets sharing the same name represent distinct binding sites for 
the given transcription factor or miR. Members of c5 gene sets share 
a common gene ontology annotation, and c5 gene sets are named for 
that GO annotation. 

Analysis of sperm miR environment 

To examine the role of paternal stress exposure in regulation of 
sperm miR environment, we assayed changes in sperm miR expres¬ 
sion following breeding of stressed and control sires. Hierarchical 
clustering analysis using Pearson’s correlation as a metric success¬ 
fully segregated pubertal and adult stress samples from control sam¬ 
ples, as visualized in a heat map (Fig. 4A). Expression of the following 
miRs was significantly increased in both pubertal stress and adult 
stress samples: miR-193-5p (F (2j9 ) = 7.3, p = 0.01), miR-204 
(F( 2 i 9 ) = 5.7, p = 0.02), miR-29c '(F (28) = 5.3, p = 0.03), miR-30a 
(Fp’s) = 17.2, p = 0.0008), miR-30c (F (2>8) = 16.6, p = 0.001), 
miR-32 (F (2>6) = 5.3 ,p = 0.046), miR-375 (F (28) = 13.5 ,p = 0.003), 
miR-532-3p (F (2 7) = 5.4, p = 0.04), and miR-698 (F (2 8) = 6.3, 
p = 0.02) (Fig. 4B). Genes predicted to be targets of at least three of 
the nine significantly increased miRs are shown in Table 4. 


Discussion 

Dysregulation of stress neurocircuitry is a common underlying fea¬ 
ture in neuropsychiatric disease (Arborelius et al., 1999; Moghad- 
dam, 2002; Nestler et al., 2002; Walker et al., 2008; Corbett et al., 
2009) . Adverse parental life experiences likely contribute to offspring 
HPA stress axis dysregulation through germ cell epigenetic modifi¬ 
cations (Bale et al., 2010; Franklin et al., 2012). Previous studies have 
demonstrated that transgenerational metabolic, endocrine, and 
stress phenotypes can transmit through the male lineage following 
an initial maternal exposure, findings that emphasize the dynamic 
reprogramming of male germ cells (Jimenez-Chillaron et al., 2009; 
Franklin et al., 2010,2011; Guerrero-Bosagna et al., 2010; Dunn and 
Bale, 2011; Morgan and Bale, 2011). Furthermore, identification of 
offspring physiological and behavioral changes following direct pa¬ 
ternal perturbations has highlighted sperm epigenetic marks as a 
potential mechanism of transmission (Carone et al., 2010; Ng et al., 
2010; Dietz et al., 2011; Vassoler et al., 2013). Therefore, in our study, 
male mice received 6 weeks of stress before breeding, and offspring 
outcomes were evaluated. As epidemiological studies suggest that 
the greatest effects on offspring and grand-offspring health occur 
during developmental windows of increased germ cell vulnerability, 
we exposed males to stress either throughout puberty or only in 
adulthood (Kaati et al., 2002,2007; Pembrey et al., 2006). To identify 
potential epigenetic marks in sperm associated with the transmission 
of offspring phenotype, sperm miR content was assessed. 

In our HPA stress axis assessment, we found that offspring from 
both paternal stress groups displayed significantly blunted cortico¬ 
sterone responses to acute restraint. While maternal stress has been 
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Table 4. Common predicted gene targets of the nine identified sperm miRNAs 
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Gray boxes indicate that the given gene was a predicted miR target in miRanda, miRDB, miRwalk, and TargetScan 
alorigthms. Asterisks denote a predicted function in chromatin regulation, DNA methylation, or miR processing. 


shown to elicit sex-specific responses in stress-related outcomes, pa¬ 
ternal stress surprisingly affected stress reactivity similarly in both 
sexes, suggesting distinct modes of maternal versus paternal trans¬ 
mission (Mueller and Bale, 2007, 2008; Franklin et al., 2010, 2011; 
Kapoor and Matthews, 2011; Morgan and Bale, 2011). It was also 
unexpected that HPA stress axis dysregulation occurred regardless of 
the timing of paternal stress exposure. Though epidemiological 


studies of Swedish famine years highlight preadolescence as a period 
of heightened germ cell susceptibility to malnutrition, our data illus¬ 
trate that germ cell reprogramming by stress exposure occurs in 
adulthood (Kaati et al., 2002, 2007). Relative germ cell vulnerability 
throughout the lifespan may depend on the type of perturbation 
encountered, with severe caloric restriction and chronic stress elicit¬ 
ing distinct mechanisms. 

The specificity of HPA stress axis reprogramming by paternal 
stress was remarkable, as is highlighted by the lack of significant 
differences in any of the assessed behavioral measures. These 
outcomes support divergent functioning and regulation of physio¬ 
logical and behavioral stress reactivity. Stress pathway dysregulation, 
either increased or decreased reactivity, reflects an organism’s inabil¬ 
ity to respond appropriately to a changing environment, and thus 
may characterize a disease state or a predisposition toward one. Al¬ 
ternatively, a reduced physiological stress response may reflect an 
adaptive response programmed by the paternal lineage as a protec¬ 
tive measure, ensuring greater offspring fitness in what is expected to 
be a more stressful environment. Interestingly, rodents models have 
previously demonstrated offspring HPA axis dysregulation follow¬ 
ing maternal or paternal exposures to diverse challenges including 
alcohol, lead, chronic stress, or social defeat, suggesting that the HPA 
axis is especially sensitive to environmental manipulations—the 
proverbial canary in the coal mine—thus warranting greater mech¬ 
anistic examination (Coe et al., 2003; Weaver et al., 2004; Virgolini et 
al., 2006; Hellemans et al., 2008; Kapoor and Matthews, 2008; Mu¬ 
eller and Bale, 2008; Dietz et al., 2011; Harris and Seckl, 2011; Mor¬ 
gan and Bale, 2011). 

To assess the neural circuitry underlying the blunted HPA in 
paternal stress offspring, we first used a pharmacological approach of 
an acute administration of the SSRI, citalopram, to stimulate the 
HPA in the absence of external stress. The serotonin system is a 
potent modulator of HPA stress axis reactivity, and has been re¬ 
ported previously to be disrupted in offspring following maternal 
stress (Heisler et al., 2007; Mueller and Bale, 2008; Franklin et al., 
2011; Huang et al., 2012). However, we found that citalopram 
evoked similar corticosterone production and recovery levels in con¬ 
trol and paternally stressed offspring, suggesting that the blunted 
response detected following a restraint stress did not likely involve 
dysregulation of serotonergic projections to the PVN. Given intact 
serotonergic regulation following paternal stress, principal compo¬ 
nents of the HPA axis, the pituitary and adrenal glands, were also 
examined in paternally stressed offspring for altered gene expression 
that may explain stress pathway dysregulation. We found no expres¬ 
sion differences between groups in the pituitary gland (CRF 
receptor-1, the CRF receptor on pituitary corticotropes, and pro¬ 
opiomelanocortin, the precursor to ACTH) or in the adrenal gland 
(melanocortin receptor-2, the ACTH receptor in the adrenal cortex, 
and 11|3HSD-1, the rate limiting enzyme in corticosterone produc¬ 
tion), suggesting that the reprogramming of the offspring HPA axis 
by paternal stress exposure lies centrally. 

By microarray analyses, we examined broad changes in gene ex¬ 
pression in key stress regulatory brain regions, the PVN and BNST. 
We hypothesized that paternal stress would more likely affect tran¬ 
scriptional patterns, rather than single candidate genes, as a complex 
cascade of programming events separates an initial paternal pertur¬ 
bation from its eventual effects in the adult offspring brain. Gene set 
enrichment analyses provided a characterization of these patterns by 
analyzing expression of a priori defined gene sets based on known 
biological associations, e.g., shared biological pathway, conserved 
cis- regulatory element, or common gene ontology category 
(Mootha et al., 2003; Subramanian et al., 2005). The enrichment of 
gene sets related to HPA stress axis regulation provided intriguing 
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potential mechanisms that may underlie shifts in offspring stress 
responsivity. Most notably, enriched PVN expression of glucocorti¬ 
coid receptor responsive genes could confer tighter control of CRF 
neurons in the PVN in paternal stressed offspring. Additionally, the 
majority of enriched gene sets we found in the PVN and BNST were 
those where member genes shared transcription factor or miR reg¬ 
ulation. These broad differences in gene expression suggest potential 
involvement of broad chromatin changes and epigenetic program¬ 
ming. In particular, enriched BNST expression of gene sets associ¬ 
ated with the activity of CREB, histone acetyltransferases such as 
CREB-binding protein, and six members of the miR-154 family, 
which are transcribed from a maternally imprinted locus, support 
the hypothesis that epigenetic mechanisms are likely involved in 
reprogramming of the offspring EIPA axis circuitry (Seitz et al., 2003, 
2004). 

Paternal models of transgenerational epigenetics offer the advan¬ 
tage of examining germ cell transmission in the absence of con¬ 
founding maternal factors, which can make deciphering the timing 
and identification of parental epigenetic transmission difficult. 
Among the limited set of epigenetic marks carried in sperm, miRs 
are poised to be products of a dynamic paternal environment and to 
uniquely contribute to postfertilization gene expression. Therefore, 
we assessed the sperm miR content from paternal controls and stress 
groups. Elierarchical clustering of over 250 miR transcripts segre¬ 
gated all stressed sires’ miR expression profiles into a single clade, 
supporting recent evidence that sperm miRs are amenable to envi¬ 
ronmental perturbations (Li et al., 2012; Marczylo et al., 2012). Fur¬ 
thermore, the finding that sperm miR content exhibited similar 
changes in both stress groups was consistent with the effects of pa¬ 
ternal stress exposure on offspring stress responsivity, and therefore 
supports a continuous germ cell vulnerability to epigenetic repro¬ 
gramming across the lifespan. 

In examination of significant changes in specific sperm miRs, 
paternal stress both throughout puberty and in adulthood signif¬ 
icantly increased the expression of nine (miR-193-5p, miR-204, 
miR-29c, miR-30a, miR-30c, miR-32, miR-375, miR-532-3p, 
and miR-698). As miRs are normally present at low basal levels in 
sperm, the observed variability in control miR expression and the 
directional increase in miR levels with stress were not surprising. 
Importantly, miRs have been shown to regulate gene expression 
after fertilization, and may impact offspring development by act¬ 
ing on oocyte stores of maternal mRNA (Giraldez et al., 2006; 
Rassoulzadegan et al., 2006; Grandjean et al., 2009). Assessment 
of the top predicted mRNA targets of the nine significantly in¬ 
creased miRs revealed DNMT3a (DNA methyltansferase 3a), a 
critical regulator of de novo DNA methylation important for im¬ 
printed genes, and two proteins involved in miR processing, 
trinucleotide repeat containing 6b (Tnrc6b) and metadherin 
(Mtdh) (Kaneda et al., 2004; Lazzaretti et al., 2009; Yoo et al., 
2011). These genes, targeted by four of the nine upregulated 
miRs, suggest that paternal stress may affect offspring HPA axis 
by intervening in the epigenetic regulation of oocyte develop¬ 
ment. Alternatively, sperm miRs may have altered epigenetic 
programming earlier, during spermatogenesis (Amanai et al., 
2006). 

In summary, our findings demonstrate that male exposure to 
stress, either throughout puberty or in adulthood, reprograms pa¬ 
ternal germ cells and results in transmission of an offspring HPA 
stress axis dysregulation phenotype. Whether such diminished stress 
reactivity would be detrimental or beneficial to offspring likely de¬ 
pends on the environment into which they were born, as well as 
genetic background factors. However, the finding that mild stress 
experience can produce such long-term changes in male germ cells 


provides an important and novel mechanism contributing to neu¬ 
ropsychiatric disease risk. 
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